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Abstract 

Here I present a new discrete model of quantum mechanics for rela- 
tivistic 1-electron systems, in which particle movement is described by a 
directed space-time graph with attached 4-spinors, but without any con- 
tinuous wave functions. These graphs only consist of few space-like edges, 
e.g. the ground state of atoms is described by two nodes and one edge, 
and interactions only take place at the nodes. 

The fundament is an extremal principle for a relativistic invariant 
"LAGRANGian sum" , from which "field-equations" and "equations of mo- 
tion" are derived, so the states (including the graph nodes) are completely 
determined. 

As important validations of the model, the corresponding graphs for 
the stationary DlRAC-equation for the atom are drawn and the correct 
spectra are computed (Sommerfeld- levels). 

Also a discrete SCHRODINGER approximation and an associated 
"HAMlLTONian sum" are derived and the correct equation of a classical 
moving particle under LORENTZ-force is presented. 

I hope, that this new approach will help, to overcome some problems 
of current quantum mechanics by making the wave function superfluous. 
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1 Introduction 

In this paper, I propose a new discrete view to the quantum world, without the 
use of a wave function concept. 

The wave function was introduced by Erwin Schrodinger around 1925 
to describe quantum mechanical states, like electrons inside an atom, for which 
classical descriptions failed. 

However, there is a long, unceasing discussion about the interpretation of 
this wave function, especially for the measurement process ("collapse" of wave 
function) (see e.g. JU| or pp. 40, for a comprehensive discussion). 
On the other hand, it seems paradox, that the description of discrete quantum 
states (like energy levels of an atom) required the invention of a new continuous 
field. This new field has the additional strangeness, to be 'not physical', i.e. is 
not directly measurable, like all other known fields. 

Additionally, the current QED-theory has severe difficulties arising from 
infinite integrals, which have to be eliminated by some mathematical tricks 
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(renormalization). Many physicists believe, that these are at least suspicious 
(see e.g. [8], p. 458, 0, pp. 166). As far as I see, the theory presented here, 
will not show any infinite values. 

As strong motivation to try discrete theories, 1 want to cite A. Einstein 
from one of his last works (see [5], Appendix II, 1954, p. 163): "Man kann 
gute Argumente dafiir anfuhren, daB die Realitat iiberhaupt nicht durch ein 
kontinuierliches Feld dargestellt werden konne. Aus den Quantenphanomenen 
scheint namlich hervorzugehen, daB ein endliches System von endlicher Energie 
durch eine endliche Zahl von Zahlen (Quanten-Zahlen) vollstdndig beschrieben 
werden kann. Dies scheint zu einer Kontinuums-Theorie nicht zu passen und 
muB zu einem Versuch flihren, die Realitat durch eine rein algebraische Theorie 
zu beschreiben. Niemand sieht aber, wie die Basis einer solchen Theorie 
gewonnen werden konnte." 1 2 

Nowadays, there exist several proposals to introduce 'discreteness' into 
physics. Most of them postulate a space-time lattice at the PLANCK-scale, of 
about 10 _33 cm and 10 _44 s. These scales are assumed to play a fundamental 
role in general relativistic quantum gravitation (which is not considered in this 
article). However, due to the smallness of these units, it is not to expect to 
find consequences of the lattice structure with currently available measurement 
techniques. The quantum fields in these theories mostly appear as continuous 
approximations of discrete lattice fields. 

My approach is different to the above mentioned, since it considers 
discreteness at particle wavelength scales (~ 1/m for time-like edges, 
COMPTON- wavelength) , 3 i.e. much larger, so it is directly related to the 
quantum nature of the particle. 

On the other hand, I do not describe the whole space-time as a gridded 
structure, only the movement of the elementary particles should be considered 
as not being continuous but in finite steps. There might exist an underlying 
PLANCK-scale grid, but this is not needed in the following considerations. 

In this paper, I deal with a new view to special relativistic quantum 
mechanics of spin-1/2 particles in electromagnetic fields, i.e. Dirac equation 
and their solutions (in flat Minkowski space-time), where the wave function is 
a 4-dimensional complex field. 

In the discrete theory proposed here, the moving particle 4 is described 

1 Emphasations by A. Einstein 

2 Translation: "There are good arguments, that reality cannot be represented by a contin- 
uous field. It seems to follow from quantum phenomena, that a finite system of finite energy 
can be described completely by a finite number of numbers (quantum numbers). This seems 
not to fit to a continuum-theory and must lead to the attempt to describe reality by a pure 
algebraic theory. Nobody sees, however, how the basis for such a theory can be achieved. 

3 In this article I use "natural units" , where h = 1 and c = 1 holds. 

4 Of course, resting particles can be described as special stationary cases. 
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as a set of space-time points, with finite space extent and, of course, infinite 
time extent. Attached to each time-like edge is a DiRAC-spinor (as a constant) 
which replaces the continuous spinor field. 

It turns out, that the number of points required to model e.g. the energy 
levels of the atom, is only in the order of the quantum numbers. In this (sta- 
tionary) case the time-edges of the graph are simply equal and the space edges 
are constant. In nonstationary cases, however, also graphs with bifurcations 
and combinations are imaginable, but these are not considered in this article. 

The equations for both grid-points and spinors are derived from an extremal 
principle of one general sum-function. This seems to be the most appealing 
aspect of this new theory. All other discrete theories, known to me, postulate 
some preset, fixed grids, which arc not influenced by the fields. 5 
This extremal principle resembles a Lagrange functional, which is widely 
used in quantum physics and esp. quantum field theory. Here, the LAGRANGian 
is replaced, of course, by a sum over the space-time graph and the variational 
principle simply maps to the variations of the points and spinors. 

To show the correctness of the theory two important cases are discussed 
and computed below: the stationary electron in the CouLOMB-field (atom) 
and an electromagnetic acting particle in nonstationary case (accelerated by 
LoRENTZ-force). 6 

A last word to the structure of this article. Some of the evaluations are 
not strictly needed in this paper, esp. the simple start cases. However, it is to 
expect, that many readers are not very familiar with the unusual notations used 
here. Thus, I think it is always better to start with the simplest possible cases 
and then proceed to the more complex states. 

In any case, I tried to put the evaluations straightforwardly as possible. So 
many proofs are left out, or shifted to the numerous footnotes, so that quick 
readers can skip them. Longer ones were put into the appendix, which thus 
became quite voluminous. It also serves, to illustrate the correspondencies to 
classical theories and formalisms. 

2 Notations, Entities and Transformations 

In contrast to the usual description of DiRAC-spinors as 4-spinors, I use a slightly 
different notation by 'spinor-matrices'. These are complex 2x2-matrices, i.e. 
they have the same number of components. 7 

°In some vague sense, this resembles the concept of general relativity, where the space-time 
metric g^ v is influenced by the mass distribution. 

6 Since most formulas derived there, are simple algebraic, they can easily be implemented in 
numerical computer programs. In fact, I have done this for most of the examples, to validate 
the evaluations numerically. 

7 They may be thought of writing the two bi-spinors as a two-column matrix, see Appendix. 
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MiNKOWSKi-vectors and LORENTZ-transformations are then represented as cer- 
tain subsets of these 2x2-matrices, with constraints explained below. 

In the appendix it is shown, that both notations are equivalent for the 
DiRAC-equation and classical relativistic electromagnetism. 

The main reason for using this form is, that all entities are represented by 
the same algebraic structure, and many of the following equations are much 
better readable, than in component notation. 

General 2x2-matrices are here denoted with uppercase letters P,Q,S, 

The usual operations with the matrix P = ( a '\) , with complex a, b, c, d, here are 
written as: 8 

• p = (_ d c '^) : adjuncted matrix of P, 9 

• \P\ = ad — be : the scalar determinant, with PP = I |P|, 

• T(P) = a + d: the scalar trace, with P + P = I T(P), 

• P T — : the transposed, P* = ( a ,' d ») the complex conjugated, 

• pt = (P*) T = : the adjungated (or HERMiTEan conjugated). 10 

Since it is often needed in the following, and not quite obvious, I state 
here the general circularity relation for the trace of any matrices ABC, . . . , X: 
T(A BC ...X) = T(BC ...X A). 11 

A Minkowski- vector is in this formalism represented by a HERMiTEan ma- 
trix, and here denoted by boldface (upper- and lowercase) letters: 12 = M, 
to distinguish it from other matrices (spinors, transformations, electromagnetic 
field tensor). 

It has, of course, 4 real components, which can be mapped to space-time coor- 
dinates (i, x, y, z) in the following way (see e.g. [H], pp. 16): 13 

ft + z,x ~ iy\ , . 

M= ' y )=tl + xa 1 +ya 2 +za 3 , (1) 

\x + ty, t — z) 

where Oi are the usual PAULi-matrices. 



8 / denotes the 2x2- identity matrix. 

9 The inverse matrix of P is then, of course P _1 = P/\P\. 

10 A11 operations here commute, e.g. (Pt) = (Py and products obey {PQy = Q* P^ and 
[PQ) = QP. 

n It can be derived from the symmetry relation T{AB) = T(BA), which again follows e.g. 
from \A + B\ = {A + B)(A + B) = \A\ + \B\ + T(AB) = \B + A\. 

12 One exception is the relativistic 9-operator. 

13 Since every matrix M can be uniquely decomposed into a HERMiTEan and anti-HERMlTEan 
part by M = A + iB, with At = A, = B (which then transform independently un- 
der LORENTZ-transformations), all 2x2 matrices can be seen as generalization of MlNKOWSKl- 
matrices. 
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lorentz- transformations are represented by unimodular matrices 
T, \T\ — 1, and thus have 6 real degrees of freedom. 14 Ordinary space 
rotations additionally fulfill the condition = T, leaving 3 free real degrees, 15 
while special LORENTZ-transformations obey = T. 16 

With respect to their behaviour under space-time transformations, we must 
distinguish between spinor- and MiNKOWSKi-matrices. 

Let T be a LORENTZ-transformation, then a spinor transforms with P — > TP 
(then follows e.g. pt — > P t T t ), while a Minkowski matrix transforms with 
M -> TMT'. 17 

The determinant |M| is then obviously the MiNKOWKSKian invariant 16 (al- 
ways real): 

|M| = t 2 - x 2 - y 2 - z 2 . (2) 

It is remarkable in this formula, that the signature of the metric tensor (H ) 

automatically follows from the property of HERMiTEcity. 

To build the general scalar product of two Minkowski matrices A, B serves 
the formula, which is obviously also invariant and real: 

T(AB) = T(BA). (3) 

From the representation QJ should also be noted, that the trace of a 
Minkowski matrix maps to the time component, and the operation of "ad- 
junction" is a space (P 3 ) inversion. 



3 Space-Time Graph and "LAGRANGE-Sum" 

3.1 Descriptions of Particle Movement in SRT 

This short section is intended to explain concisely, how particle movement is 
described in the context of Special Relativity and Minkowski space-time. Also 
a discrete variant of movement (which is not used in the following) is sketched, 
but quantum effects are not considered here. 

A continuous relativistic particle trajectory is given by a space-time curve i.e. 
the 4 functions x(r) = (x(r), y(r), z(r), t(r)) 19 which is usually parametrized 

14 Mathematically speaking, in terms of LlE-group theory, T build the SL(2, C)-group, which 
is a double cover of the LORENTZ-group. 

15 By this definition they build a subgroup (the quaternion group, see appendix), while 
special LORENTZ-transformations do not. 

16 E.g. a matrix T = (ol'^a) w ith real f3, performs a (t, z) transformation. 

17 The above condition for space rotations =T = T -1 then leads to M — » TMT" 1 , 
consequently the trace of T(M) = 2t is invariant, as required. 

Consider e.g. the transformation T = iai, which performs a rotation of 180? around the 
x-axis. A full rotation is then represented by T = (io"i) 2 = — 1. 

18 The proof of invariance is simple: |M| -> |T||M||Tt| = |M|, since \T\ = 1 

19 Or in matrix notation written as: x(r) = f'^l+^f [T^7 iyl ; T I) ■ 

V ' \x(T) + iy{T),t(T) — z(T)/ 
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by the eigentime t. 20 

This t is defined by dr 2 = \dx\ = dt 2 - dx 2 — dy 2 - dz 2 . 

Since dr 2 is an invariant, all Minkowski- vectors dx can be classified by its 
sign: dr 2 > : time-like, dr 2 = : light-like, dr 2 < : space-like. 

Usually by the condition of causality it is required, that no interactions over 
space-like separated regions occur. The movement of a particle is restricted to 
time-like vectors (resp. light-like for massless particles). 

The discrete form of such a space-time curve is then simply a sequence 

xi = (x!,yi,zi,ti), x 2 = (x 2 , V2,z 2 ,t 2 ), ... 

This sequence can be considered as a graph with the edges (xfc — > Xfe + i), 
which describe a movement in finite "jumps" and "time-likeness" here obviously 
means |xfc + i — Xfc| > 0. 

The (continuous) movement of a classical charged particle in an electromag- 
netic field (LORENTZ-force) can also be derived from a variation principle for a 
LAGRANGian . 

Let all possible space-time curves be parametrized by a parameter A : x(A). 
In matrix notation the action integral is then written: 21 



^2 

£(x(A)) = J 



dx _. - . . dx. 

Tx +eT(A(x W 



Ai 



The variation of space-time curve x(A) = x e (A) + Sx(X) results in an extremal 
curve x e (A). 

def 

Then A is identified with the eigentime of the extremal curve dr = ||dx e ||. 

The extremal curve is then given by the equation of LORENTZ-force: 22 

d 2 x e ( dx „ x dx\ 
m ^ 2 -2W + F Tr)- 

It is interesting to state here, that there exists also discrete representations of 
the above LAGRANGian formalism. However, this example is given as illustration 
only, and not used in the following sections! 

A possible discrete variant of the above integral is: 

C = ||xfc+i - x fe || +e^T(A(x fe )(x fe+ i - x fe _i)). 

k k 



20 But also other parameters, e.g. t may be used. The advantage of using r is, that the 
elocity vector u = ^ is then normalized to unity. 

21 absolutc value of a MINKOWSKI vector written as ||x|| d = \/jx[ 

22 See Appendix for the relations between A and F in matrix notation. 
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Herein the variation of one space-time point (node) x^ leads to a discrete version 
of the LORENTZ-Force: 

/ x fe+ i-x fc x fc -x fc _i \ ef + \ 

m T n - rr ~ o ( x fe+! - x k-i)F k + Fl (x fe+ i - Xfe_i) . 

V||x fe+1 -x fc || ||xfe-x fc _i||y 2 V V 

Like in the following sections, here also the idcntifaction ||xfc+i — Xfc|| = 1/m is 
possible, since this expression is conserved (approximately). 23 

3.2 General Considerations with Space-Time Graphs 

Space-time grids are commonly used to solve partial differential equations, e.g. 
numerically. Then integrals (e.g. the LAGRANGian functional) are represented 
as sums. Usually, the gridded structure is viewed as approximation of the 
continuum, and the smaller the edges are, the better the approximation. 
In this theory, I try another point of view: the grid represents the quantum 
state and the differential form is the approximation. 

In fact, it turns out in the following, that e.g. to describe bound states in the 
atom, that there exist "minimal grids", which suffice to represent the exact 
states in the DlRAC-theory 

It should be noted, however, that the usage of finite (esp. space-like) edges 
introduces some kind of nonlocality and causality violation into the theory. 24 

I will start with the general expression for the "LAGRANGian sum" over any 
graph. Let {x^} be the nodes as MINKOWSKI matrices (numbered in an arbitrary 
order), Hy,Mj be some MiNKOWSKian matrices, A the electromagnetic vector 

def 

potential, which has the values Aj = A(xj) at the grid points, and Pi some 
spinors: 25 

C = ^T((x, - x^Hy) + e£T(A 4 M i ) - 2m£>(|P i |). (4) 

ij i i 

However, the first double sum is not to be applied for all pairs (i, j), but only 
for edges. 

The auxiliary matrices , M j in this equation shall be constructed as 
HERMiTEan bilinear forms from the fundamental spinors Pi, as explained below. 

At first, however, it is to prove, that this sum fulfills all requirements for a 
"LAGRANGian" : it is a real scalar and invariant under all lorentz transfor- 

23 This has no analogon in the continuous case. To prove it, multiply the eq. with q = 
Xfc_|_i — xj c _ 1 (e.g. from left) and take the trace (i.e. building scalar product with q), Then 
the r.h.s. vanishes, since 7(F) = 7(F^) = 0. 

Also should be added, that the used graphs show some similarities to FEYNMAN-graphs. 
However, bifurcations are not considered here, and the mathematical background is completely 
different. 

25 As usual, m denoting particle mass, e electrical charge and SftQ the real part of a complex 
number. 
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mations. 

It is obviously scalar, by construction. To prove the reality of C, I state that: 

• for any HERMiTEan matrix A holds trivially: T(A) = real, 

• for any two HERMiTEan matrices A, B holds: T(AB) = real, due to the 
symmetry relations: T(AB) = T(BA) = T(B+At) = T((AB)+). 

Since all factors are HERMiTEan matrices and 3?(|P|) is always real, the 
complete sum is real. 

To prove lorentz invariance, I state, that the epressions x _1 and A 
transform with T^(-)T and therefore the expressions T(x _1 H) and T(AM) arc 
invariant scalar products. 

The determinant \P\ is trivially invariant under LORENTZ-transformations, if 
the spinor transformation rule P — > TP is considered, q.e.d. 



3.3 Regular Space-Time Graphs 

Now I consider regular space-time graphs. The restriction to these graphs is 
mainly due to the problem, that it is not yet clear, what physical conditions 
can lead to bifurcations or combinations, and the mathematical difficulties in 
handling them. 

This is no principal limitation, and as shown in the following, many problems 
of one-particle quantum mechanics can be described with these graphs. 

One first introduces double indices for the nodes {x^}, where the first index 
should stand for space, while the second index k stands for time steps (thus 
unbounded, k = — oo • • • oo). 26 

The regularity condition then means, for any time index k there exist 
n nodes: i = l,...,n and that for any k should hold |xjfc — xjk\ < 
(space-like edge), and for any i,k: (xj^+i — Xj^l > (time-like edge) and also 
ti.k+i — U,k = ^1~{ x i.k+i — x i,fc) > (direction of time arrow). Also, only 
timely consecutive nodes (k, k + 1) shall be connected by an edge. 

The (constant) spinors Pi are now considered to be assigned uniquely to the 
time-like edges (not to the nodes): (x^+i, x^) <-> Pa-. 

Of course, this assumption introduces a fundamental asymmetry between space 

and time and leads to different formulas for the above introduced H. 

For space-symmetry reasons, the following ansatz is suggested: 

H i;fe +i,i,fe = PikP} k , H i;fejife = Pj,fc_i o P-j.k, 27 (all other combinations of 

26 This numbering scheme does not violate the LORENTZ-covariance of the following evalua- 
tions. 

27 The circle stands for the HERMiTEan conjugated expression: P o Q d = i(PQt Qpt). 
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indices have no edge assigned) and M.ik = ^{PikP} k + Pik-iP^-i)- 28 
Then from (J3J results: 

c = j2 T ^ + i-^M)- i p lk pi) 

k.i 

+ e^T(i(A iik+1 + A i;k )Pi fc 4) - 2m^5R(|P lfc |). 



(5) 



A: , i 



To visualize the kinematic terms (first and second term) of this sum, the 
following picture is used, where the spatial extent number is set to n = 2 (this 
example graph e.g. also represents the ground state of an electron in an atom): 
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Fig. 1: example space time graph (time axis vertical) 

The extremal principle now considers the sum £ as a function of all inner 
variables Xik^Pik, whereas possibly some boundary variables have to be fixed, 
to account for initial conditions: 



£(-x. ik ,Pu 



Extr. 



(6) 



3.4 Simplest Case 

To demonstrate the method of deriving the "field equations" and "equations of 
motion" from this principle, I start with the simplest case: no electromagnetic 
potential (A = 0) and the graph has only one spatial index: n = 1. 
This model represents a freely moving spin-1/2 particle. 

The graph then reduces to a sequence of MINKOWSKI space-time points, 
which is actually a discrete particle trajectory, as explained in the section T3. II 

28 all obviously HERMlTEan 
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The spatial index can be omitted, and the second and third sum in equation 
JSJ) are zero. It remains the sum: 

C(x k ,P k ) = ^T((x fe+1 -^ k )- 1 P k Pl)-2mY,^(\Pk\)- (7) 

k k 

This sum is obviously invariant under the "local" transformations P k — > P k S k , 
when S k Sl = \S k \ = l. 29 That means, the spinors P k are determined only up 
to these factors by the following equations (gauge invariance). 

At first, I consider the variation of one specific P k , in this sum. This variation 
is similar to the usual methods in quantum theories. 

For the variation of the real part of the determinant, is used: 29?(|P|) = |P| + |Pt| 
and 



|P + 5P\ = (P + SP){P + SP) re |P| + SPP + PSP = \P\ + T{6PP), (8) 
therefore results from substitution P k — > P k + SP k the variation 

SC. = T((x fe+1 - Xk )-\6P k Pl + P k 5Pl) - mT(SP k P k + 8P£a%). (9) 

To simplify the formulas, I define an auxiliary variable d = (x^+i — x*.) -1 , so 
equation writes (after using the general circularity relations for the trace): 

5C = T{5P k {Plv k - mP fe )) + T(6Pl(v k P k - mpj)). (10) 

As usual, the variations of 6P k and 6P^ are considered as independent, therefore 
both terms must vanish. An expression T(XY), however, can only vanish for 
any matrix X, if Y = holds. So the two (equivalent, since by definition v = vT) 
equations result: 

PfeV fe = mP k and w k P k = mP\. (11) 

This equation corresponds to the usual "field equations" in quantum mechanics. 
Here it forces (by taking the determinant on both sides) that, since |vfe| = real, 
also \P k \ = \p£\ = real, consequently |vfe| = m 2 = const, or \x k+1 —x k \ = 1/m 2 , 
implying that the motion vector is a time-like vector of the constant length 
1/m. 30 

The second variation 31 (which has no correspondence in current theories) 
varies the nodes x k . For this, it is needed to state, that for small & << x holds 
(see appendix) 32 

{x + Sx.)- 1 wx" 1 -x _1 &x _1 . (12) 

29 These S are normalized quaternions, see appendix. 

30 In the rest frame of the particle, thus trivially holds At = 1/m, for all others relativistically 
At > 1/m. 

31 This is not to misinterprete as a second order variation. 

32 A simple proof is, to multiply the equation from the left (or right) with (x + &). 
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The variation of influences only and v k ~i (by their definition) and leads 
to 

<5vfc = +Vfe& fc v A . and dvk-i = -v k _ 1 &c k v k _ 1 (13) 
The variation of C in is consequently (second line by circulation) : 

SC = T{v k dn k v k P k Pl) - Tfvfc-i&fcVfc-iflfc-xP^i) 

- T{6x k (v k P k Plv k - Vfc-xflk-xP^jVfc-i)). (14) 

Again, this expression must vanish for arbitrary leading to 

VfePfcP^Vfc = vnPn^tjVn, (15) 

Inserting equations (|llfl twice, results in -Pfc-P/t = Pk-iPl-^ = const.. This is 
fulfilled by the condition P k — P k —%S, with an arbitrary matrix, that obeys 
SS^ = I (gauge invariance). 

On the other hand, follows from equation (fT5)) : 33 

Vfc=v fc _i, (16) 
that says, that the motion vector is constant over time, as it should be. 

The result of above computations is, that the particle "moves" in jumps 
with a constant time-space vector Xfc + i — x k = Ax, which is related to the 
particle wavelength by yj\ Ax| = 1/m. The Minkowski vector v% — const, is 
therefore to identify with the relativistic energy- impuls vector p = e + p, where 
e denotes the energy and it holds |p| = e 2 — \p\ 2 = m? . 

The spinor orientation does not have any influence on the motion, as to expect 
in the absence of an external field. 

In any case, the time-steps that arise in these jumps are by orders too small 
to be visible in experiments (e.g. with ultrashort laser pulses). For an electron 
e.g. holds At w l(T 20 s. 

3.5 Stationary Case 

This case describes bound states, e.g. an electron in an atom. I will show in 
the following, that it leads to the correct energy spectrum. 

The grid for this case is considered as "time invariant" , i.e. it is claimed 
Xj,fc+i — Xjfc = t, where r = real, t > is a constant time step. 34 
Only in this "periodic" case, the space components of the graph are constant 
over time (for all time indecees). 

33 The other solution v fc = — v fc _ 1 would imply a step backwards in time. 
34 Likc above, it can be identified with the inverse energy of the state and its value follows 
from the equations below as the eigenvalue. 



A new discrete view to quantum mechanics 



13 



The space-like edges shall be pure space vectors (traceless matrices x = — x). 

Additionally, for all i, k the ansatz Pn-+i = Pi,kS is made 35 (i.e. 
P%,k = Pi.iS^ 1 ) with a constant matrix S, obeying SS^ = I and \S\ = l. 36 

It follows |P i)fc | = |P 1)fc |, P hk+1 P[ k+1 = Pi, k Pl k = ■■■ = /', ,/',;, = const. 

and r,j, ,/•;, = r,j, t s'r; h . ••• Pi,iSfiplJ 

The summands of C in (JSJ then become independent of the time index if 
also the external field is considered as time invariant: A^/. = Af i. The index 
k = 1 can be dropped, also the summation over k can be omitted and one gets: 



- 2m^»(|P i |). (17) 

i 

Again, the double sum is to build only over edge-pairs 

To simplify the formulas, I introduce a set of auxiliary variables 

def , _ _ _ 

Ujj — (X^ Xj J ^ 

The antisymmetry of the factor in the double sum in i,j leads to a simpli- 
fication, e.g. the summands for the pair (1,2) are: 

Py&P} + P 2 SPl - P 2 &Pl - PiSPl = P^ - S)Pl + P 2 (S - S^)P}. 

Any unit quaternion S can generally be represented with real A and U as general, 
pure vectorial, unit quaternion (W — U — —U, \U\ — l) 37 as: 

S = e xu = cos A + U sin A. (18) 

This gives S — = 2[/sinA. Inserting this in (|17l) follows, that the double sum 
is proportional to sin A (which is the only term containing A). The extremal 

principle then requires (since A is a free ansatz-parameter) , that cos A = 0, 
therefore S — ±U and S* = —S, 38 and equation (|17fl simplifies to (here also r 
is replaced by the energy e = 1 /r) 

C = E n(e + ek^pp}) + E nvuPjSP}) - 2m £ (19) 



35 Considering S beeing independent of the spatial index i is the standard method of "sepa- 
rating variables" (here time and space are to separate). In function form one would make e.g. 
the ansatz ty(x,t) = (f>(x)ip(t) and here S k stands for the time dependency factor. 

36 These conditions hold for all unit quaternions S, which is equivalent to S £ SU(2). 

37 U has 2 free real parameters and its general form is U = ( 1 sm > cos v e \ 

Vcos tpe 2 sm Kp/ 

38 In the following one may use e.g. S = U = (n'— i) i the explicit form does not matter. 
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The variation of Pi is carried out like in the last section and leads to the "field 
equation" (a system of n linear equations for the Pi, i — 1, . . . , n): 39 

(e + eAi)Pi + J2 uuPjS = mP} (20) 

ij 

This equation becomes equal to the DiRAC-equation for the stationary case, 
if the "operator" u ij is replaced by the spatial derivative (V-operator). 
This correspondence is shown in the appendix. 40 

Again, I consider a second variation of the grid points x,. The "equations 
of motion" derived with this method, result in the determination of the grid 
points. Again, this procedure has no counterpart in the present theories. 
Please consider again, that this variation does not affect the spinors Pi . 
In the next section I will show, that it produces the correct quantum states for 
the electron in CouLOMB-potential. 

For the variation of jq in the sum (|19fl . at first it is needed, that 41 

<5A 4 = <5A(x 4 ) = -TC&iflJA,, (21) 

where d is the MiNKOWSKian differential operator (an explicit representation in 
matrix notation is given in the appendix). 
The variation of Uy = (xj — Xj)^ 1 is again 

Siiij = —Uij&CiUij and Siiji = ■ u,,i i k,u J ,. (22) 

It results for the variation of x,-, when another auxiliary variable Hij = PiSPj 
(with h\a — —Hji) is introduced for simplification (note, that in the sum i is 
fixed) 42 

SC = ^T(,k 4 a)r(A l P l F/) + ^ 7 iu.^w,,//,, - u, /! k,u,,//„ :. (23) 

ij 

From the demand 5C = for all &i (with operations like in previous section), 
results the equation 43 



e 



dTiAiPiP}) + ]T ujiiHij + //; n h , = 0. (24) 



39 Note, that the sum in this equation, in contrast to the sums above, is simple, since i is 
fixed. 

40 If one consideres this equation as a classical eigenvalue-problem for e (given Xj and A;) 
it has, of course, at least n solutions Ej and corresponding eigenvectors (at most 4n, because 
every matrix-equation has actually 4 scalar equations). 

41 see appendix for a short explanation. 

42 the differential operator here operates only on the external field A, of course, since Pj is 
considered as a constant. 

43 for the important special case, that A is time constant, (d + 8) A = holds and therefore 
from 1241 can be derived „ u^T^Hij + Hj^) = 0, since u| 4 = scalar. 
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The combined solution of (|2U[) and l|24[l is then the expected quantum state, 
which also determines e, as it is to see in the following examples. 

However, sometimes it is easier to use another method, that obviously leads 
to the same results. For all solutions of equation (|20|l follows, that £ = holds. 
Then C(e,p\,p2, . . .) = is an implicit equation for e {pi subsumming all free 
variables, here as simple reals). Since C shall be extremal with respect to all 

de 

other parameters pi it follows — — = 0, so e(»i,...) itself must be extremal 

dpi 

(usually minimal). 



4 Electron in an Atom 

In the usual approximation, the atom nucleus shows a COULOMB potential, lead- 
ing to scalar eA = eV(r) — 2f, where a — e 2 w 1/137 denotes the "finestructure 
constant" and r = \J x 2 + y 2 + z 2 the EUCLiDian distance to the center. 44 

4.1 CoULOMB-potential 

To solve the equation l|2(J|) we first rewrite it for COULOMB-potential, by using 
the auxiliary parameters £j = e + a/rj, and the same with conjugation 45 

e l P l +Y^ "ij Pj S = mP} equals EiP} - ^ u i3 P] S = mPi . 
By dchning 

Pt= f Pi+^ and Pr*f (Pi-pl)S (25) 
and addition/subtraction of both equations one gets: 46 

(e l ~m)P+ = -J2mjPf and ( £i + m)Pr = £ Uij -P/. (26) 
The LAGRAGIAN sum <(T§|) reads with this substitution: 47 

c = E m )i^ + i - & + m )i p ri) + E r (W^ ■■ (27) 

i ij 

At this stage the SCHRODiNGER-approximation is easily feasible, by setting 
£j+msi 2m in l|26|) or (|27|l . This is done in sectional However, in the fol- 
lowing sections I want to present precise results. 

44 In matrix notation holds r = ^—\x\. 
45 consider u = — u, ut = u and St = S. 

46 Since P+ and P~ by their definition obey: (P+)t = P+ and (P~)t = -P~ this is a 
quaternionic decomposition of P as P = P + — P~ S, where P + and iP~ are quaternions (see 
appendix). Since also iu is a quaternion, this set of equations 1261 can represented with these 
algebraic entities. 

47 Of course, one also can derive again 1261 from 1271 . 
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4.2 Ground State 

At first, I start with the simplest case: the ground state in an atom, and will 
show, that the above formulas lead to the correct energy (and atom radius). 
This state will be described with only two points (n = 2) 48 : xi,X2. Then 

exists only one edge U12 = (xi — x^ -1 == u = — U21. However, as explained in 
the section IF1 this edge is counting twice, giving a factor of 2, so equation Q26JI 



reads: 49 



(ei - m)P 1 + = -2uP 2 ~ and + m ) p i = +2uP, 



+ 
2 ■ 



(e 2 - m)P+ = +2uPf and (e 2 + m)P 2 ~ = -2uP+. (28) 

In these four equations only the two spinor-pairs P x + , P 2 and P^ 1 ", Pj - are cou- 
pled, therefore eliminating the P, - , results in the pair of equations: 

(ei - m)(e 2 + m)P+ = 4u 2 P+, (e a + m)(e 2 - m)P 2 + = 4u 2 P+. (29) 

In these two equations the spinors P^~ are then freely variable and can be 
divided out (except for the two singular cases P± = P 2 ~ = or P 2 + = P^~ = 
resp.) 50 and one gets 51 



(ei — m)(e 2 + m) — 4u , {e\ + m)(e 2 — m) — 4u . (30) 
Consequently follows from (ei — m)(e 2 + m) = (ej + m)(e 2 — m), that ei = £ 2 

,1 ii- ^e/ 

must hold, i.e. r± = r 2 = r. 
The resulting equation is 



(e + ~Y - m z = 4u z . (31) 

Since |xi| = |x 2 | = — r 2 holds, by triangle formulas one gets: — |xi — x 2 | = 
4(r 2 — h 2 ) < 4r 2 , with h < r as height on the edge, so 4u 2 = l/(r 2 — h 2 ) and it 
results 52 

<r,h) = -^± ] fJ7^^. (32) 

As explained above, the LAGRANGian extremal principle requires that this ex- 
pression is to make stationary with respect to r and h, which after simple com- 
putations gives 53 immediately h = and r = ^ 1 ^ a2 1=3 (bohrs formula for 

£_ ° ^_ am am v 

48 It is easy to show, that for n = 1 no stationary solution exists. 

49 with Ti d = ||Xj|| = j^ii and again e t d = e + f- 

50 The above decomposition on P into P ± in @U has the consequence that the LAGRANGian 
1271 becomes a sum of two terms: C\ (P^~ , , e, xi , X2) + £2(^2^, P\ , e, xi, X2), where the 
spinor-pairs must be viewed as varying independently. Therefore, these singular cases not 
have to be considered. 

If one nevertheless computes these cases in full detail, eg. the first one, where only the second 
cq. of 1301 holds, it turns out, however, that the resulting expression e(xi,X2) has only one 
stationary point, which is no extremum, but a saddle-point. 

51 The factors on both sides are simple reals. Consider again, that u is a traceless matrix 
by definition, so u = — u and u 2 = — |u| > 0. 

52 In this formula r, h should not be misunderstood as usual variables: they get fixed values 
after using the extremal principle, also e is then a constant, of course. 

53 only the + sign of the root gives for positive a (attractive potential) an extremum 
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the radius of the hydrogen atom) and finally the correct energy of the ground 
state: 

e = my/ 1 - a 2 . (33) 
4.3 Space-Grid for the General State 

The purpose of this section is, to present a general space grid, that is stationary 
together with the spinors. And I will show, that this represents the correct 
quantum states of an atom. 

Therefore, I first consider a separation of variables, namely the radial 
variable r and angular variables (on the sphere). This separation is possible due 
to the symmetry of the COULOMB potential and is similar to usual procedures. 

However, I want to emphasize here, that it should be possible to find 
more general methods, which do not rely on the assumption of separability, 
used below. I think, that once the edge scheme is fixed (i.e. which nodes are 
connected by edges), it can be proved, that the nodes are general stationary 
points. This was e.g. shown for the ground state in the last section, where it 
was forced by the equations, that both nodes have the same distance to the 
center. 

The grid should consist of n spheres with the radii r; = r\,..,r n , and all 
these spheres should have the same set of node normals. That means that every 
point of the grid can be represented as = TiPj, where Pj are unit vectors: 

|Pil = -l- 

Two points on different spheres / Tj also should only be connected by an 
edge, if they have the same spherical coordinates. 

Again, the task is to find a stationary point for the LAGRANGian sum in (|27(1 : 

C(n, Ph Pt,Pr)- 

The separation ansatz now assumes, that also the spinors Pf 1 can be factorized, 
namely as 

P+ = ftAi and Pr = g iPi Ai, (34) 

where /$ and gi should be real constants only depending on the radial index 
and Ai matrices, only depending on the angular index. 

Considering that on radial edges holds Xj — Xj = Pi (r{ — rj), introducing the 

auxiliary parameters ef 1 d = €{ ± m and inserting the above into the sum l|27|l 
leads to (the summation labels R, S should denote summation over radial, resp. 
angular indicees): 

VfceS / \i£R ij£R % 3 J 
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J2 T(A k p k ( Pk - Pi )"%) ) ^ (35 ) 

The usual separation idea is now, that the angular and radial dependent factors 
in both summands must be separable, that requires with some constant n: 

2k 53 1^1 = £ TiAiPiipi-p^Aj). (36) 
s s 

Then the sum can be decomposed into two independent factors C = CsCr, 
where 2k arises as eigenvalue in Cs and the radial factor Cr becomes: 

c* = E( e r/f + 4 9 f) - 2*x; s -f- + e f ': Jj - w 

' i ... 'i '7 

The solution of the angular part Cs is given in the appendix, section IF1 

4.4 Solution of the Radial Equations 

By variating the Ut9% m one gets: 54 

a^ft-Er 1 - and ^i = -/i+Er^-- (38) 

3+1- 3^1 

At this point, I want to emphasize the correspondence to the radial dif- 
ferential equations, derived from DiRACs-equation, with similar presumptions, 
namely they read: e~ f = ~g - g' ', e+g = ff + f (see e.g. 0). 
The detailled discussion of this and also the connection of the both associated 
LAGRANGians is given in the appendix. 



55 



The second variation, which considers the ri, additionally gives (for all ri) 

The equations (|38|l set up a system of linear equations, which can be 
considered as eigenvalue problem for e (if all are fixed). Together with l|39|l 
they form a set of 3n equations for the 3n + 1 variables r^, fi, gi, e. 56 
Since the first system is linear in fi,gi and the second bilinear, however, they 
are normalizable, consequently the number of equations equals the number of 
variables, indicating that only discrete solutions exist. 



54 Note, that the sums are only over the radial index, from here on. 

55 Consider, that de^/dri = —a/rf and the double sum contains each term twice. 

5B From here on, n denotes the number of spheres, not nodes. 
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For n = 1 the solution can be derived directly, giving (correctly) 57 

e = — \/ n 2 — a 2 , n = \/ k 2 — a 2 . (40) 



am 



For the general case, it turns out, that a simple linear ansatz for the fi,gi, 
where a, 6, c, d are real constants 

fi=a + bri and gi^c + dn, (41) 
and a set of fj, obeying the equations (with two parameters A,7) 58 

V^ = A-2 (42) 

Z — ' r i - r 3 n 

gives a solution which results (after longish computations) in the correct formula 
for the energy levels, where n r = n — 1 is the radial quantum number (see e.g. 
0, p. 126): 59 



(?) 



2 



a 



- / = 1 + ( ai> ¥— ■ ( 43 ) 

5 Electron under LORENTZ-Force 

This section is intended, to demonstrate the working of the method for one 
non-stationary case. 60 

I consider here the same case as in eq. (0, except that also an electro- 
magnetic field is present, i.e. again with n = 1 (no space-like edges) and 

v fe = (x fc+ i - x fc ) L -y L 



C = ]TT((v fc + e -{k k+l + A k ))P k Pl) - 2m£R(|A|). (44) 

k k 

The "field equations" become similarly: 

(v fe + ^(A k+1 + A k ))P k = mpj. (45) 



57 Consider from |3SJ: e 1 fi = -frgi, and e+gi = ^-/i and from 1^1 ^(/j 2 + gf) = 2/t/igi. 
This system of 3 equations has only the above solution. 

58 These equations are related to LAGUERRE-polynoms, and discussed in the appendix. 

59 The parameters are then determined as 7 = \/ tt 2 — a' 2 and A = V m' 2 — e' 2 . 

60 To tackle problems of this type as initial value problem (e.g. numerically), one should 
consider the following method. Suppose the three consecutive points xi,X2,X3 are used 
(and the associated spinors P\,Pz). Then the variation equations for the inner point X2 and 
P\iPli 1451 1471 are used, building a set of 3 implicit equations. In this set, the initial values 
of xi, X2, P\ are inserted, which eventually results in the values for X3, P2, and so forth. 

ei Again this sum is obviously invariant under the "local" transformations P k — * P k Sk, when 
S k Sl = 1 and \S k \ = 1. I.e. the spinors P k are determined only up to these factors by the 
following equations. This is a partial analogy to gauge invariance of standard DIRAC theory, 
except the vector field A is not transformed here. 
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By taking the determinant on both sides of this equation immediately follows, 
that \P k \ must be real 62 and consequently it must hold for all k (it is closely 
related to the conservation of energy) : 

|v fc + ^(A fc+ i + A fe )| = m 2 = const. (46) 

One remarkable consequence of this simple formula is, that regardless of the 
history of a particle, in case of vanishing vector potential, in its rest frame 
always holds At = 1/ra. 63 

The second variations (for Xfc) result in: 64 

v fe P fc P,!v fe - Vk-iPk-iPL^k-i + - A dT{A k (P k Pl + Pk-iPL,)) = 0. (47) 

Please, consider again, that the set of equations l|45|l and 147f) must be solved 
simultaneously. 

By multiplying l|45(l from the right with Pi and (v + eA) -1 from left one gets 65 

P k Pl = rnlPfcKvfe + ^Afc+i+Afc))- 1 

« m|P fc |(v- 1 -v- 1 ^(A fc+1 +A fe KT 1 ). (48) 

Inserting this and the corresponding term for P k -\P}._ 1 in (|47|l and omitting 
terms ~ e 2 A 2 gives (after division by m): 

|P fe |(v fc - |(A fc+ i + A fc )) - \P k -x\(v k -i -|(A fc +Afe_i)) 

+ ^T(A fc (|P fe |v fe - 1 + |P fe _ 1 |v^ 1 ))=0. (49) 

Reordering gives (after a bar-operation): 

|P fc |(v fc -|(A fe+1 + A fe ) + I^A^ 1 )) (50) 

= \P k -i\(vk-i - |(A fc + A fe _ x ) - |9T(A feVfc -^)). 

Now with 1 = Xfc +1 — Xfc (by definition) the approximations are used: 66 

Afc +1 « Afc + ^T(v- 1 B)A k and A k _, » A fc - ^(v^^Afc (51) 



62 consider |v + eA||P| = m 2 \P\* and |u| = real for any HERMlTEan u. 

63 The equation, however, reveals some important new issues: 
E.g. adding a constant offset to A (which does not affect the classical LORENTZ-force) here 
changes the discretization and thus modifies the results. It seems, that apparently no full 
gauge invariance can be derived for this model. For small fields, however, the results are equal 
to the classical theory. 

Then for example, consider the simplest case: a resting particle in scalar potential A = 
t/(x). Then eq. gHJ reads: + U) 2 = m 2 i.e. At = Since At > is supposed, 

only the range — oo < U < \m\ for the external field is possible. 



64 consider again (SAj. = iT(<5xfcc9) A 



65 if the approximation of small field |eA| << |v| is used. 
66 see appendix, chapter "Differential Calculus" for explanation. 
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resulting in 

|P fe |(v fe -eA fc t T {y-^B)A k ) + |9T(A fe v^)) (52) 

= |i\_i|(v fc _ 1 -eA fc + P -T(v-\5)A k e -dT(A k v-\)). 

Now the equation for the heldtensor F (at point x^) is used (see appendix |B|) . 
which for any u obeys: uF + F^u = T(du)A — dT(Au) giving 

|P fe |(v fe - e A, - ^(v-^fe + ^v- 1 )) (53) 
H^-iKvfc-i-eA* + l(^-i F k+FH-i))- 

Since one cannot generally claim |Pfe| = |Pfc-i| (which case could be easily 
solved) , the symmetrical ansatz (which is always possible, of course) with a new 
real variable A: \P k | = 1 + A, |Pfc-i| = 1 — A is used. Also, the centered difference 

q d = v^ 1 + v^j = Xfc + i — Xk-i is used. Then one gets: 67 

(1 + A)v fc - (1 - A)v fc _! = ^(qF fc + i^q). (54) 

To determine A, this equation is multiplied (from right) with q = v^ 1 + vl_, 
and then is taken the trace, so the right-hand side vanishes (since F + F = 0). 
It remains 

T((l + A)(l + vfev-^) - (1 - A)(l + Viv, 1 )) = 0. (55) 

Then, with a d = — v^^i and another auxiliary matrix B d = VfevZ"\ = 
1 + av^j, one gets: 

(1 + A)r(l + B) = (l-A)r(l + J B- 1 ) i.e. X = J+t^ItIbY (56) 

With T{B- r ) = and \B\ = jr^ 68 and the approximation T(B) w 2 

follows A 

To state the approximation of LORENTZ-force of eq. I|54|l it remains to ex- 
plain, that the relativistic velocity vector u = dx/dr = x is discretized as 69 
u = (xfc+i - x fe _i)/^|x fc+ i - x fc _i| « ^q i.e. q w and a is actually a 
discretized acceleration vector: 



a = v fc -v fe _i = |vfc|(xk+i - x fc ) - |vfc_i|(xfe - x fc _i) 

^ x fc+ i - 2x fc + x fc -i _ cPx 
Z_ (At) 2 ~ dr 2 ' 



(57) 



6 since |A| <C 1 is supposed to be small, only the dominant term ~ A is considered, which 
is v 

6s consider from eq. K + f (A fc+1 + A fc )| = [v fc _i + f (A fc + A fc _i)| (= m?) and 

|eA|< |v| thus |v fc | « |vfc_i| holds. 

69 centered around x fc , using At = y/\x k+1 - XfcZiJ = ^/iv^ 1 + v~\\ ss yj 4|v~ x | S3 2/m. 
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Finally from eq. (|54|l results the equation of LORENTZ-force (see appendix) with 
small corrections: 70 

a=^-(uF fc + ^u)-2-u. (58) 

2m m 



6 SCHRODINGER-Approximation and "HAMILTON- 
Sum" 

This section shall demonstrate, that also a discrete form of the classical 
stationary SCHRODINGER equation and its associated HAMiLTONian can be 
derived as approximation from the above discrete dirac formalism. This 
approximation is always possible for electrons in weak electromagnetic fields. 71 

I start with equation l(2l)|) from section ETT1 which represents a bound state 
of an electron, but here in a general electric potential field V(x) (real scalar, 



time independent) with Vi = V(xi) 

{e + Vi- m)P+ = -J2 ».../', and (e + V t + m)Pr = u.^ + - (59) 

3 3 

As mentioned, herein the SCHRODiNGER-approximation is easily feasible, by 
setting e + Vi + mtn 2m in the second equation. 73 Then P~ can be expressed 
directly with it: 

^«iE u ^ + - (80) 



.72 



2m 

3 



Inserting this in the first of H59J) . gives with E d = e — m as energy: 74 

(E + V t )P+ =-^-J2 u Mk P k- (61) 



2m 

jk 



Please note, that the double-sum is to build over all edge-pairs and (j, k). 
This is still a matrix equation. To get a scalar equation from it, one adds the 
adjuncted, resulting in a scalar expression: 75 



(E + V,){Pt + Pt) = -7^ »M',//V + P^kjujil 

jk 



'"The correction term resembles a corresponding term in DlRACs motion equation, which 

2 __ 
reads in this notation and scaling: 6 e 7rm (a — (see e.g. p. 173) 

71 The way of deriving this approximation from DlRAC equation is similar to standard QM. 

72 and again uy = (x; — x_, ) _1 for the space-edges 

73 This is the usual approximation method for small energy, since e ~ m and V <C m 
74 in many textbooks the potential energy U = — V is used instead of V in the formulas 
75 consider that E, Vi are scalars and Uj^ = — Ujk and Ukj = — Ujfc holds. 
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and then drops the vector part of P^ in the r.h.s. (ignoring all spin-effects), by 
setting it a real scalar ipk = P k = P k - 

The result is a discretized form of the stationary SCHRODiNGER-equation: 

(E + Vi)i/)i = ^2 -(uijUjk + u kj Uji)ip k . (62) 

jk 

The sum ^2~2jk( u ij u jk + u kjUji)ipk there represents the second order partial 
derivation operator Aipk 76 at the point x«. 

One can also easily define a " HAMiLTONian-sum" , from which the above 
SCHRODiNGER-equation (Mi2t can be derived again (by variation of all ipk)'- 77 

H{i>i, x*) = J2( E + + UiiUjfc^V'i (63) 

i ij k 

This sum obviously corresponds to the classical HAMiLTONian for the stationary 
case. 78 

6.1 General Considerations about Ground States 

At first, I want to discuss the solution for the ground states of any potential V, 
which is described here with n = 2 nodes. 

Then one has only one edge with u d = u i2 = (xj — x 2 ) _1 = — u 2 i and two 
equations: 

1 u 2 u 2 

(E + Vi)^i = -— Ui 2 u 2 i^i = — i>\ and (E + V 2 )ip2 = -^—^2- (64) 
2m 2m 2m 

2 I 

It follows E + Vi = E + V2 = j—, i.e. at first one can conclude V\ =V-2- 
Please note, that for every MiNKOWSKian pure space-vector u (with u = — u) 
holds u 2 = — |u| > 0, therefore E > -V1.2 

Also notable is the decoupling of tf)\ , ip2 m equation (|64|) , meaning that their 
values are independent. 79 

However, in contrast to the case of CoULOMB-potential discussed below: if 
there exists a stationary point xo of the potential (with dV g* ^ — 0), then also 
a solution with only one node (n = 1, no edge) is possible, which then would 
have the energy E = — V(xq). 

These solutions do not have counterparts in current QM. To exclude them, 
there must be a principle, that excludes stationary solutions with n = 1. At the 
moment I cannot see, what this can be. 

7e It is always scalar and real, of course. 

77 Of course, it it also possible to derive this HAMiLTONian directly from the LAGRANGian 
sum, with the same assumptions. 

78 The proof, that the triple sum is always a real scalar, is simple when again the hermitecity 
and antisymmetry of the Uij = ut^. = — uji is used, e.g. (ui2U23)^ = U23U12 = U32U21. 

79 This does not hold for the exact solution, given in chapter 4.2., there the spinor P k cannot 
vanish at any point k. 
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6.2 Ground State in CoULOMB-potential 

The CouLOMB-field (of an atom nucleus) is (with r = | |x| | as euclidian distance) 



n 



V(x) = -. (65) 
r 

Then (again from V\ = V2) directly follows r\ = r2 = r. As explained above, in 
contrast to the one-dimensional case, the edge is counting twice again, so one 
gets 

^ ot 4u 2 

E=— + — 66 

r 2m 

By setting 1/u 2 = (xi - x 2 ) 2 = 4(r 2 - h 2 ) results: 80 

The condition §§ = then gives h = 0, and §f = gives r = 1/am (atom 
radius) and finally the energy of the ground state of hydrogen: 

£7= -my. (68) 
6.3 Quantum Harmonic Oscillator 

This section is included, to give readers the most simple testcase of the discrete 
theory. The one-dimensional harmonic oscillator has the field (with usual scale 
factor): 

V (x) = -^lo 2 x 2 . (69) 

Here for simplification is set m = uj = 1, giving from l|63(l the "HAMiLTONian 
sum": 81 

H( fc «)-B*-i^4s; (70) 

ijk 



Here all node-pairs shall be connected by one edge, with a total of 
edges. 

The "field equation" (by variation of t/j m ) results in n eqn. 

2{E r 2 ^ m + E fx -x% - Xi ) = °- (71) 

^ V m ^31x^3 ^kj 

3K 

The variation (differentiation) of (|15() by x m ,m — 1, .., n gives also n eqn: 

80 With simple triangle formula for the triangle (xi, 0, X2), with /i as height and h < r. 

81 Consider e.g. the one-dimensional case with x = (^'o) = 2 = Then the Ujj = 
(Xj — Xj) -1 all commute, and their products UijUj^ = ^ x ._ x . jr^ -x k ) are a l wa y s simple 
scalars 
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9 



2-m VVn 



2 ^ (x m Xj) 2 {xj Xk) 2 ^ (x.i x m )^(x m ^/c) 
1 v - Vi^fc 1 v - VjV"™ 



2 ^ (x.j .T m ) (x ni Xk ) 2 _ (x.j ) x m ) 2 
This can be simplified to: 82 

- z ro v4 - E ^ + E = ( 72 ) 

These 2n equations (|71l I72|l are to solve with 2n variables (ipi,Xi). 
Since it is wellknown, that HERMiTE-polynoms are the eigenfunctions of the 
classical quantum harmonic oscillator, it is suggested using them to find a so- 
lution. And indeed, it is simple to prove the following solution with their help. 83 

I will show in the following, that with Xi as zeros of these Hermite- 
polynoms, and the most simple ansatz for the ipi~. ipi = a (an arbitrary 
constant) all 2n equations are fulfilled, so this is actually a stationary point of 



set 84 



The zeros of HERMiTE-polynoms (x%, . . . ,x n ) obey the implicit equation 

" 1 

E — = x * ( 73 ) 

This is easy to prove with the methods given in appendix E: "Orthogonal Poly- 



The proof of l|72l) is then simple (factors a 2 dropped) with reordering both 
double sums: 



El / X - J \ -> I \ _ \ -> X m — Xj 
fW T-l^ I L. 4 m ™ . L. 4 ™ qn. / Z 4 { W T'\^ 



(x m Xj) 2 y ^ X m Xk ^ Xj Xk J ^ (x m Xj) 2 

From the above implicit sum formulas (|73|> for the Xi, one can easily derive 86 

n n 

} Xi = and > — = x 2 . — (n — 1). 

Xk — Xi 

i—l i^tk 



82 by the reassignement of the sum indicees in the 4. sum i —> k (equals — 1. sum) and of 
the 3.rd sum i <-> k (equals — 2.) and then i — > j 

83 Possibly this is not the only solution, also if permutations of the Xi are considered. 

84 Please note, that the X{ are uniquely determined by 1731 (up to permutations, of course). 
The first few zeros can be computed explicitely, e.g. for n = 1 : x\ = 0, for n = 2 : X\,2 = 
±y/~^, for n = 3 : xi = — ^J\, X2 = 0, x$ = ^f\- 

85 They have the generating differential equation (see eq. I A. 301 1: y" — 2xy' + Xy = 0, i.e. 
u(x) = 1, v(x) = -2x, Z = -2x. 

86 e. g . with y. = y . = yv, t-i + -^-) = -< n - n + x\. 
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The double sum in l|71|l then becomes: 

{x m — Xj)(xj — Xk) (x m — Xj) 

and one gets: 

E=^ (74) 



This formula reproduces the wellknown energy levels of the quantum harmonic 
oscillator, if the restriction n = even is made. 87 

However, it is not clear yet, which physical principle excludes the other solutions 
(for n = odd). 



7 Conclusions and Outlook 

Here I presented a new discrete view to quantum mechanics, where the 
continuous wave functions are replaced by a space-time graph with attached 
constant spinors and the differential equations by discrete, algebraic equations. 
These equations are derived from a general "LAGRANGian sum" over the graph. 

The remarkable new idea is, that the graph nodes are not to be arbitrary 
set, but determined by the variation principle for the same sum. 
Since the graph is a MINKOWSKI space-time graph, this includes the time steps 
(which are for stationary cases then determined by the energy eigenvalue of the 
state) and gives also a valid description of particle movement (nonstationary 
case). 

With this model many classical problems of quantum mechanics are solved 
and give the expected results. (However, some of the solutions do not have a 
correspondence in classical QM. It is not yet clear, which physical principle 
suppresses them.) 

It is thus my hope, that this model can make the wave function obsolete 
(similar to the light aether, that became obsolete by the theory of Special 
Relativity), and all quantum phenomena can be described by a finite number of 
numbers, as an algebraic theory, like A. Einstein suggested. 

However, there remain also many unanswered questions: 
One principal task is, to introduce real dynamic behaviour into the theory, e.g. 
to describe emission and absorption processes. Then also the model of a photon 
should arise. The space-time graph for a photon, as massless particle must be 
described differently, however, since the time-like edges must be replaced by 
light-like edges, with |Ax| = 0. 

These processes could be probably modelled with bifurcations and combi- 
nations of the graph. In general, due to the implicit character of the formulas, 

87 Consider the units % = u) = 1 
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this should be possible, because they may have more than one solution. 

Additionally, the usage of the electromagnetic vector potential in the theory, 
can only be seen as an approximation of interactions with (virtual) photons. 
This however, would imply significant changes of the variational principle. 88 

Another interesting aspect is the question, if it is possible, to set one 
primary entity (between spinors and MiNKOWSKian vectors), from which the 
other can be constructed as derivation. Apparently, this can be only the spinor 
part. 

However, it is to expect that this question can only be solved in a more general 
framework, which I suppose to be a discrete theory at the PLANK-scale level 
(quantum gravity), that will have a quite different concept of space-time. From 
that, the presented theory will arise as approximation. 

Other important tasks are: 

• Is it possible, to simplify the LAGRANGian sum 10, e.g. combine time-like 
and space-like terms, and find some kind of deeper explanation for it? 

• Can gauge invariance represented better? What happens for strong fields, 
where | |eA| | ~ m? 

• How are many-particle systems described? This should be possible, of 
course, by defining a composed LAGRANGian . 

• How are antiparticles described in this theory? 

• Is it possible, to embed this theory into the framework of general relativ- 
ity? 

At last, it is to ask, of course, how far the theory is consistent with the 
current experimental knowledge. Especially, the representation of particle waves 
(de Broglie- waves) and entangled quantum states would be a challenge. 

A DlRAC-Equation in Matrix-Notation versus 
usual Spinor-Notation 

As stated above, I will show here, that both notations are equivalent. 

For this purpose, I start with the conventional representation for 4-spinors: 89 

ij^ip = mip. (A.l) 

88 However, it should be stressed, that for the most important case (couLOMB-potential), 
there exists a striking correnspondence between the factors of the kinematic terms in sum CT . 
i.e. (xi — Xj) _1 and the potential term eA; = ar7 (with ri = ||x^ — xo||), that suggests, 
both terms may have the same basic cause. 

89 see e.g. HJ, pp. 24 
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There ip is a 4-column vector, and the 7 M are 4x4 matrices. Here, I use the 
WEYL-representation for the 7 M : 90 

Then the 4-spinor is decomposable into two 2-spinors ^ = (J) , which trans- 
form independently, but different (see below) under LORENTZ-transformations, 
and (|A.1(1 decomposes into a coupled system: 

i(-d° + a k d k )<f> = and - o k d k )^ = m$. (A.3) 

Now I define a (HERMiTEan) differential operator (a 2 x 2 matrix) 91 

9 =/ d° - fj fc a fe (A.4) 

explicitly: 

9 =(t lt^a + S V ) = ll + V - (A - 5) 

V Ox ' dy ' Ot i " 9z / 

Then the equations (|A.3|I read: 

- id$ = and - id* = m<i>. (A. 6) 

The LORENTZ-transformation T, \T\ — 1 here operates as follows on the entities: 

d -> TdT\ 8 -> f t<9T , $ -> Tt$, * _> r*. (A.7) 

Then both equations (|A. 6|) are obviously covariant under this transformation. 
As usual, an electromagnetic interaction is introduced by the substitution 
0^ — ► — ieA^, which gives here: 

(i<9 + eA)$ = -ra* and (id + eA)* = -m$. (A. 8) 

Then the second equation of (|A.8|I is converted in the following manner. 

One states the general formula for every 2x2 matrix M (M T denoting transposed 

matrix): M = crM T a, with a = ia 2 = (_°io)- 92 Since $EM>) has the form 
Al^f = — m<&, with M — id + eA, one gets aM T a^ — — m<I>, which can be 
rewritten to 93 M J a^> — — mcr$. 

Of this one takes the complex conjugate, where (M T )* = M>: 

AfW* = -mcr$*, with Aft = -i9 + eA. (A.9) 

90 The 3 PAULl-matrices are again cri = U'J), o"2 = (°'~ ! )i CT 3 = (o'-i) 
91 note, that I2 = I2 and ctj. = — <tj. holds, i.e. 9 = 9° + cr^d k and all I^^k are HERMiTEan 
matrices 

92 The reason is, that the bar-operation means spatial inversion (x,y,z) — » (—a, — y, —z), 
and that is equal to the combined operation of transposing (i.e. y —* —y) and a rotation 
around y of 180?, given by T = i<T2. 

93 trivially, since a = — a, a 2 = — 1 
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One then defines a new operator for 2-spinors \fr = a^* , which obeys 4> = — ^> 
(since a 2 — —1) and with that eq. I|A.9|I then writes (—id + eA)^ = — to$. 
One then can combine both equations into one 2x2 matrix equation 

eA($, $) + id($, -*) = -m(#, (A.10) 

Now one defines the "spinor-matrix" P d = ($, $) and states Pt = -($, $), 94 
and with the auxiliary matrix 5 = (g'_°) finally gets: 95 

eAP + dPS = mPt , (A.ll) 

Note, that according to above definitions P transforms consistently with 
P — > T P under LORENTZ-transformations and the equation (| A. 1 1|> is obvi- 
ously covariant. 



Gauge covariance is in this notation represented with the local transfor- 
mation P — > PU, with 17 = e AS (where (A(x) is a real, scalar function of 
space-time). Then holds (the parentheses are set here, to denote the action 
of the differential operator d): d{PU) = (dP)U + (d\)PUS. Then, equation 
IjA.llfl is covariant, if the simultaneous transformation eA — * eA + d A is used. 

The stationary case is given with the ansatz P(r,t) = Po(r)U(t), with 
U(t) d = e- £ts , which (since d t P = -ePS) 96 results in: 

(eA + e)P + VP S = mP f . (A.12) 



B Relativistic Electrodynamics in Matrix- 
Notation 

I will shortly sketch here the basic equations of relativistic electrodynamics in 
matrix-notation without explicit proofs. Each equation can be checked, e.g. by 
converting it to usual component notation. 

The tensor of the electromagnetic field is defined from the vector potential 
by: 97 

F d = ^(9A - Ad). (A. 13) 

J" An explicit proof is » = (£), * = (£), 5 = (jg.), P = (£•_$), pt = (£;_£), 

pT=(:^;|)=- W >). 

95 By using S 2 = — 1 follows (e.g. by TAYLOR expansion): e AS = cos A + Ssin A = (g e — ♦$)• 
96 Since the scalar operator dt commutes with Pq and dtU = —eUS. Also obviously U 

commutes with S and C/t = {/. 

97 Contrary to usual notations, here differential operators like d can operate to the right, 

resp. left. In ambiguous cases, therefore the operand should be marked. 
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It is a traceless matrix (F + F = 0, by definition) and obeys the transformation 
rule F — ► T^FT^ . It can be decomposed into a HERMiTEan and anti-HERMiTEan 
part, that are the electrical and magnetical field vectors, which both are HER- 
MiTEan, traceless matrices (E^ = E, = B): 

F = E + iB and F^=E-iB. (A.14) 

Therefore, it is obvious, that both transform independently under spatial 
rotations, but are mixed under special LORENTZ-transformations. 

The MAXWELL-equations are simply (with J as current, also HERMiTEan) 98 

J = dF, (A.15) 

and the equation of continuity (follows from last eq. with F + F = 0) reads as 

T(dJ) = dl + 3d = 0. (A. 16) 

Finally, the LORENTZ-force on a particle with mass m and electrical charge e, 
that is moving with the relativistic velocity vector u = dx/dr" results in an 
acceleration vector a = du/dr: 100 

a=-l(uf + F t u), (A.17) 
2m 



At last, I have to derive an identity for the last term of above equation, valid 
for arbitrary u, which is used in the section [SJ 101 

uF + F+u = i(u(9A-A<9) + (A<9-<9A)u) 

= i(u(T(«9A)-2A9) + (2Aa-T(<9A))u)) 

= Adu - uAd = Adu + Aud - Aud ~uAd 

=o 

= T(9u)A - 9T(Au). (A.18) 

C Differential Calculus and Approximations for 
Matrices 

The formulas stated here are standard vector analysis, they are shortly listed 
here for readers, not so familiar with the notations in this paper. 



98 these are actually 8 real equations, for the real and imaginary parts! 
"t is the eigentime, given from dr = -y/jdxj 

100 The orthogonality of a, u is written as T(au) = and follows directly from IA.17I . 



101 Note, that T(<9A) = T(9A). 



A new discrete view to quantum mechanics 



31 



The total differential for any field (matrix or scalar) U is given simply by 
dU{x + dx) = U(x + dx) - U(x) = ^T(dxd)U. (A.19) 

The simple explanation is, that the scalar invariant is explicitly written: 
±T(dx<9) = dtd t + dxd x + ■■■ 

One often above used approximation is for the expression (X + where 
X, S are both matrices, with \X\ >> \6\. However, the following approximation 
holds for any algebra. One states 102 

(X + 5)- 1 = (Xil + X^S))- 1 = (l+X^S))-^- 1 w (l-X-^X- 1 (A.20) 



D Matrix-Notation and Quaternions 

Quaternions offer an elegant method for many computations, especially on the 
unit sphere and generally with space rotations. 

They are representable by the sub-algebra of matrices, obeying = Q. 103 
The general form is obviously, with arbitrary complex a, (3: 




(A.21) 



To represent a pure space vector v with quaternions, one uses Q = iv, which 
is obviously a quaternion, since v = —v. 104 

The quaternionic units (i,j,k) are consequently equal to (i<J\, i&2, 203), giving 
(an arbitrary real s can be added, since it is not changed by rotations): 

^ fs + iz,ix + y\ . . , . , 

Q = . = s + x\ + yj + zk. A.22 

\ix — y,s — izj 

The norm (here as matrix determinant) of a quaternion is always a positive 
real: \Q\ = \a\ 2 + \/3\ 2 = s 2 + x 2 + y 2 + z 2 . 



Ordinary space rotations are directly represented by normalized quaternions, 
following from = T (i.e. also TT^ = T^T =1), and consequently the whole 
apparatus of the 3-dimensional vector space can be drawn with quaternions. 

E.g. a rotation around any (space) axis A = /iii + + /U3k (with \A\ — 
Y^, Hi — 1) by an angle A, is represented by the transformation matrix T = 
e \A _ cos ^ _|_ ^gin.x^ which is obviously a normalized quaternion. 

102 The expansion (1 — <5) _1 = 1 + <5 + S 2 + • • • can be easily checked by multiplying both 
sides with (1 — 8). It converges, if linin^oo S n = 0, which is guaranteed by |5| < 1 
103 It is trivial, that any product of two quaternions is a quaternion again. 
104 MINKOWSKI matrices with time components, however, cannot be represented directly. 
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E Orthogonal Polynoms and LAGRANGE- 
Formalism 



This chapter is intended to illuminate the general correspondencies between 

eigenvalue problems (represented by LAGRANGians) and orthogonal poly- 
i 

noms. 

In fact, this relationship was the motivation to engage in the above theory. It 
shows, that many problems of mathematical physics, described by eigenvalue 
problems, can be reduced to small sets of equations for the roots of correspond- 
ing orthogonal polynoms. 

The pair of dirac radial equations for the 1-electron atom are used here, as 
especially related example, yet there exist many other applications. 
All following refers to the one-dimensional case, however. 106 For more than 
one dimension, there will probably exist similar methods. 

In the following, I will sketch some major relations for OP. Many of them 
(but not all), can be also found in standard textbooks, but derived with 
different formalisms. 



E.l Basic Formulas for OP 

Definition: an OP of degree n : P^ (x) is the (unique) polynom, associated 
with an interval [a,b] and a weight function w{x) > 0, x e [a, b], that is orthog- 
onal to all polynoms Q(x) of degree q < n, i.e.: 107 



These polynoms are unique up to a constant factor, of course. The roots, 
however, are unique. In the following we only deal with polynoms of the form 
(x — x\) ■ ■ ■ (x — x n ), i.e. the highest coefficient is unity. These are called monic 
polynoms. With this condition they are unique. 

It follows immediately, that two OP of different degrees n, m: 

p(n) p (m) 

are orthogonal (one is a polynom of lower degree). 

In the following, the superscript of P^ and the integral boundaries [a, b] 
are omitted, however, since they are considered fixed. 



Oo Orthogonal polynoms are mainly used in numerical mathematics, e.g. to compute inte- 
grals. Their usefulness in problem solving is widely unknown, however. 
106 Thus all variables in this section are simple reals. 

107 Thus, it roots can be computed e.g. with the following n equations: J dxwP^ = 
J dxwP^x = J dxwP^x 2 = ■■■ = J dxwP^x™- 1 = 0. 




(A.23) 
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Let P(x) = Y[ i x — x i) = ( x ~ x i) ■ • • (x — x n ), where all Xi are real and 

i=l 

distinct. 



Now one defines n associated partial polynoms to P, each of degree n — 1: 
-PfcW = 11 - /c = 1, • • • , n, see e.g. 0, pp. 502. 108 

Then by definition, P is also orthogonal to all Pk- 

The partial polynoms Pi are also mutual orthogonal (J wP\ Pi = 0). This is easy 
to see, if one expands e.g. 109 P\Pi = (°(x— X2) ■ ■ -)((x—xi)o- ■ •) = P(x~ X3) ■ ■ ■. 
Similarly follows J w xP\P 2 = 0. 

Additionally, one defines partial polynoms P k i of second order (and similarly 
of higher order): 

Pki= f H(x-Xi), P M =Pik, Pkk= f 0. (A.24) 

With these definitions, I state some important equations, which can be easily 
proved with standard methods: 

PiP 2 =PP 12 , Pi - P 2 = (xx - z 2 )Pi2, (A.25) 

Pi -Pi 



p' = Pi + p 2 + • • • , p[ = P12 + Pis + • • • = Yl 



Xi Xi 
ijtl 



One now considers the master integral 

l{x 1 ,---,x n ) d = I wP 2 = / w(x - x x ) ■ ■ - 2 {x - x n ) 2 >0. (A.26) 



It is very easy to show, that this integral is minimal w.r.t all xf. At first it 
is stationary, since = —2 J wPPi = 0. Secondly, it is a real minimum, since 

(here I define the n new constants pi = J wPf) 

g = 2/^ = 2 Pi >0 and Jg-= 0. (A.27) 

One can now (uniquely) expand an arbitrary polynom f{x) of degree < n— 1 
by the partial polynoms Pi, with n constants ff. 

f(x)=J2fi p i(*)- (A.28) 



108 They are proportional to the LAGRANGE-polynoms, which are defined as 

n ra 

^fc = 11 (( x ~ x i)/( x k - ^i)), i.e. P fc (x) = n k L k (x), with 7r fe = P fc (a: fe ) = fj (a; fc - Ej) 

109 The o mark stands here and sometimes in the following for omitted factors, to make the 
products more readable. Also, the subscripts 1,2 here denote arbitrary, but different indices 
from interval 1, . . . , n. 

110 It is then minimal among all monic polynoms. 
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Then for two arbitrary polynoms /, g (of degree < n — 1) and a linear function 
fi = a + frr, easily follows (with pi p(xi)): 111 

Wflfg = ^2 Pi^ifidi- (A.29) 



In fact, from this formula follows, that every polynom of degree < 2n — 1 can 
be integrated exactly by its values at the n grid points X4. 

All OP also obey a linear, second order differential equation (with A as 
eigenvalue): 

uy" + vy' + Xy = 0. (A.30) 

Here u = uq + u\x + U2X 2 is a polynom of degree < 2, which must not have 
any zeros in the interval [a, b] and v = vq + v\x is linear. 

If then the ansatz y = P(x) = (x — Xi) • • • [x — x n ) is made and A 7^ 0, at 
the zeros obviously must hold [uy" + vy'] x=Xk = 0. With y'(xk) = Pk{xk) and 
y"(xk) = 2Pfc(xfc) X^i/fe x k -x resu lts a system of equations for the zeros: 



2V^^ + 4^ = 0. (A.31) 

^ X k - Xi U(Xk) 

The proof, that the polynom P = (x — x±) ■ ■ ■ (x — x n ) then fulfills equation 
(|A.30|) is quite simple: Since uP" + vP' is a polynom of degree < n, and has 
(following above relations) zeros at all x k and therefore must be proportional 
to P. 113 

This set of equations uniquely determines the set of zeros, and also can be used 
to set up numerical methods to compute them. 

It also can be shown easily, that the weighting function w{x) is then related 
to the pair u, v by (wu)' = wv, i.e. w = —eJ ~ dx . 



E.2 Weight Factors for OP-Integrals 



As last prerequisite, a general, explicit formula for the weighting constants pi 
is needed (which I did not found in any textbooks). To get it, I start with the 

expression (again with Ui = u(xi)): 114 

A12 = / wuP 12 P' = / wu(P? - Pi) = u lPl - u 2 p 2 . (A.32) 

1 X? - Xi J 

llx This formula is widely used for numerical integrations, however the determination of the 
coefficients pi is often quite complicated. In the following, I will show a much more simple 
way to compute them, which I have not yet found in the literature. 

112 By this representation, all OP systems can be easily classified. Some important examples 
are: LEGENDRE-, TSCHEBYSCHEFF-, JACOBI-, LAGUERRE- and HERMITE-polynoms. 

113 The eigenvalue can be easily computed from the coefficients of x n : A n = — n(n— 1)«2 — nv\. 

114 which holds for any polynom u(x) of degree < 2 
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On the other hand, by partial integration, one easily shows, that A12 = 0, 
consequently with some constant k, follows: 115 

U\p\ = U2P2 = ■ ■ ■ = k = const. (A. 33) 

E.3 Solving the Radial DIRAC-Equations with OP 

With the help of above relations, the extremal principles can be investigated. 
Here I consider, what one may call "dual eigenvalue" problems, e.g. of the type 
of the pair of radial dirac equations, where there are two functions to variate 
independently 116 (the name of the variable x is changed from here on to r and 
the intervall to use, is of course, r G [0, 00]). The LAGRANGian here has the 
general form (where a, 6, c are some fixed functions of r) 



117 



/•oo 

£(/, g) = / dr(f'g - g'f + 2afg + bf + eg 2 ) -> extr. (AM) 
Jo 

For the DIRAC equation one has to use a = — -p, b — e + ^ ~ m, c = e + j + m, 
see e.g. [S]. 

Now the polynomial ansatz is made, with a common factor ip(r): f — ipF, 
g = tpG, where F 1 G shall be polynoms of degree n — 1 (it is presumed to be 
possible, here), resulting in: 118 

/>oo 

£(F,G,<p)= dr>p 2 (F , G-G , F + 2aFG + bF 2 + cG 2 )^extr. (A.35) 
Jo 

Since ra, rb, rc are linear expressions of r, one uses the weighting function 
w(r) = ip 2 /r, so the factor in the integrand becomes a polynom of degree 
2(n — 1) + 1 = 2n — 1, and the above apparatus can be applied: 



/>oo 

£= / drw[r(F'G-G , F + 2aFG + bF 2 + cG 2 )}. 
Jo 



(A.36) 



The polynoms F, G then are represented as F(r) — J2i fiPii 
G(r) — giPi, with 2n constants fi,gi- One now uses for the deriva- 
tions F',G' the formula: 

F' = y2f iPi , with f/ay-li+A. (A.37) 



115 The value of k can be computed by evaluating the integral J wuP' 2 , then follows k = 
~{ v + (2n — l)vi)T, where vo,vi are the coefficients of v = vq + v^x. 

116 Also LAGRANGians of ordinary, second order differential equations can be represented in 
this form. Consider for example the simple integral C(y) = J{y') 2 +ay 2 — > extr., which leads 
to y" = ay. The last eq. is equivalent to the first order pair y' = u,u' = ay, which in turn is 
represented by the LAGRANGian C(y, u) = J yu' — uy' — ay 2 + u 2 . 

117 With standard variational methods, using Jif'g + <//) = [/9]q° = 0, one easily proves, 
that it is equivalent to the pair of first order DGL: /' + af + eg = and g' — ag — bf = 0. 
118 The derivations of ip cancel out. 
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Inserting this all in eq. (|A.36|I . leads to a double sum: 

/>oo 

£ = / drwy^ y [r(figj - gufj + 2af l g J + bfofj + '■(/,</ , )!', I', . (A.38) 

Jo a 

If now the previously free variables {r*} are set to the zeros of the OP for 
the weighting function w — ip 2 /r, then all integrals can be computed, and the 
expression becomes a simple sum (again defining at = a(ri), ...): 

£ = Y^ P* r *(fi9i ~ 9ifi + 2aifi9i + hff + c t gf). (A.39) 

i 

It can be shown, that the OP to use here, are LAGUERRE-polynoms, i.e. 
they belong to a weight function w(r) = e~ 2A, 'r 27 ~ 1 , where A == V m 2 — e 2 and 
7 = \J k 1 — a 2 . 119 For these, one has u(r) = r, so from eq. (|A.33|I follows 
PiTi = const. = k and I finally get the formula: 

£ = HJ2 fi9j ~_f j + Y^r'n ,.[,,,, + b t f 2 + c l9 2 )). (A.40) 

'#3 » 

This is exactly (except the constant factor k) the same as l|38|l . q.e.d. 

Again, the similarity to the starting point 1A.34|) is remarkable, however, 
a general discussion of the preconditions for this, should be left to interested 
mathematicians. 



F Lagrangian for the Unit Sphere 

Here I will shortly derive the discretization scheme for the unit sphere, that 
leads to the above discussed solutions of DIRAC equation for the atom, in 
section IP 120 

The LAGRANGian for the angular part (p, arc the points on the unit sphere), 
that was derived there, is (see l|5B)l): 

C =^2"(A iPi ( Pi - P,-)"%) ~ |^|, (A.41) 

ij i 

where the Ai are quaternionic matrices (A, = Al), attached to the points and 
the double sum is to build over all edges (pj, Pj). 



119 The ansatz-factor function ip(r) is then ip = e~ Ar r 7 . 

120 Of course, also a more concise quaternionic representation of the following is feasible. 
However, since probably most readers are not very familiar with this formalism, I prefer the 
matrix notation here. 
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As shown below, the described grid, together with Ai is stationary, i.e. it 
makes the LAGRANGian , considered as function C(pi,Ai) extremal. 

Any point on the unit sphere (x, y,z),x 2 +y 2 + z 2 = l is represented with 
spherical coordinates i?, ip by the matrix 

z, x — iy\ fcosi), sin^e - ^ 
x + iy,—zj \ sin de^ , — cos r? 

Now one introduces the conformal one-to-one mapping of the sphere to the 
complex plane 121 x G C with: 122 

X = tan -e tip = - -. A.42 

A 2 l + costf v ' 



1 — I y I 2 

Then one gets with simple trigonometric identities: z — cos$ = and 

1 /lHx| 2 ,2 X * 



x + iy = sin ■&e %v — pp£p , consequently 



" 1+|X| 2 V2X,W 2 -1 

This matrix can be decomposed with the help of spinor-matrix factors 123 Q — 
Q(x), which I define here as: 

Q(X) = K) (A- 43 ) 



.x.-i 

and the constant matrix U d = (o'_i); namely as: 124 

p = QUQ- 1 . (AAA) 

With the help of this decomposition formula it is easy to express the required 
difference vector of two arbitrary points pi, P2 on the unit sphere (the inverse of 
this difference is the crucial part in computing the LAGRANGian sum in l|A.41|l ') 

and Qi= f Q(x&- 125 



Pi - P2 = QiUQp - Q 2 UQ2 L = Q2 (Q2Q1U ~ UQ2Q1) Qr (A.45) 



ief 
= D 2 



121 This map is the inverse of the usual "riemann sphere" map. Please note, that the complex 
number \ then transforms linear fractionally under ordinary space-rotations, namely if T = 
( Jg» „») is a rotation (with |T| = 1), then follows \ -p*£+ a * ■ 

122 E.g. the north pole (0, 0, 1) is mapped to the origin of the complex plane \ = 0, the south 
pole (0, 0, —1) to the infinite point \ = 00 and the equator to the circle |x| = 1. 
123 These factors transform similar to spinors under rotations, not MINKOWSKI space vectors! 
124 For an explicit proof consider \Q\ = -(1 + |x| 2 ), Q 2 = 1 + |x| 2 and Q" 1 = 1+ ^i Q 
125 Consider from above Q^ 1 UQk = QkUQ^ 1 
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With V 2 i d = Q2Q1 = (l + *l xi i'it~ji) it is to see, that the difference in the 

V X2 Xl-- t ~rX2Xi' 

brackets of this expression D21, is the anticommuting part of the factors U and 
V21 and can easily be computed as: 



D21 



With this equation (|A.45ll is easily invertable and one finally gets: 

1^/0, 



Pi(pi - P2)- 1 - QiUQi 1 Q 1 D^ 1 1 Q 2 = -QA ' X °- Xl )Q2. (AM) 

xi-xV 

Now one makes the following substitution for the A4 , to simplify the double sum 
in (|A.41(I . with the complex constants fii, Vi. 12& 



Ai=Q A^-^y (A - 47) 

Then follows \A\ = ifaj-IHHj ~ nv*) = ] ^l\^f ■ 

The LAGRANGian in (|A.41fl then simplifies to (it is obviously real as required): 

£{ ^ Xi) = ? 2K ? ^t^t- (A - 48) 

As always, one has to find stationary points in parameter space, which is here 
{lH,Vi,Xi},i= 1, 

Since all parameters are simple scalars, this can be done by setting the partial 
derivations zero, which gives in equations, k = 1, . . . n: 127 

BC _ dC _ dC __ () 

dn k dv k dxk 

Namely, (by the variations of v* , fi* ) result the first 2n equations (for k — 
l,...,n): 128 

yU^l 2K ^^ and ±-2*— (A.49) 

The variation of all %fe leads to the equations (k = 1, . . . , n): 129 

-Et^ + 2 fi ^|M.O. (A.50) 

Y (Xfc-X t ) 2 (1 + IXfel 2 ) 2 

Now I will show, that the point in parameter space, given by the simple 
formulas jttj = c, Ui = c\* with real constant c 130 and any set of points {xk} 

126 This is the most general ansatz, if the quaternionic restriction A; = A\ is considered. 
127 Us usual the complex conjugate of any parameter can be considered as independent, and 
since the expression C is real, the derivation by it leads to an equivalent equation. 

128 Consider the complex differentiation rules, e.g. = ^ = x* ■ 

129 The double sum contains each pair i, k twice. 

130 The resulting matrix is then simply Ai = Q~ cQi = cl. 
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(on the unit sphere by definition), that obeys 

\- 1 + xixi ■ 



0, k = l,...,n (A.51) 



fulfills all 3n conditions above. Like all others, the above sum is not to build 
over all i ^ k, but only over the edges {XiiXk), which are presented in the 
following. 131 This grid of points Xk can, as to expect, be derived from spherical 
harmonics (this is described below). 

Simple rearrangement gives two equivalent sets of equations, wherein I is the 
number of summands (edges i, k) 7 which should be equal for all Xk, k = 1, . . . , n: 

T =1 X } „ and Y Xi =-l 1 --r. (A.52) 

Then all in conditions ljA.491 IA.50|I arc fulfilled, i.e. the stationarity of 

£(/ii, i>i, Xi) is proven, if the eigenvalue n is set to K = |, q.e.d. 132 

A second stationary point is obviously given with the same set of Xk, but 133 

I 

Vi = cXt, »i = c, K= ~2' 

Please note, that both cases k = ±^ describe different quantum states in the 
dirac equation (see e.g. [8], pp. 119). 

Finally I will shortly sketch, how a set of {xk} and the assigned edges, that 
fulfill i|A.51|) . can be constructed using spherical harmonics. 
There might exist also other grids, which give the same result, but I was not 
able to find anyone. 

The nodes are supposed to be arranged on h latitude circles (which are 
defined here by |x| = const.), in the way that on every circle are 2m equidistant 
nodes. 

I.e. one has n = 2m x h nodes and they can be assigned as (with real tk > 0): 
X = t k e m ™, k = l,...h, j = l,...2m. (A.53) 

Edges shall be only on latitude circles (denoted as i £ B(k)) and longitude 
circles (as i G A(k)): 

T,^—= E ~^~.+ E — ^- (A.54) 

i Xk Xl ieB(k) Xk x% ieA(k) Xk Xl 



131 These equations are indeed covariant under space rotations, described here as linear frac- 
tional transformation of all Xi> like mentioned above. 

The simplest possible example set is given with n = 2, i.e. two points xi> X2 with one edge 

between them. Both 1A.511 . for k = 1,2, then simply require 1 + X1X2 = 0. This condition 
says, that the two points must be antipodes on the sphere, while one of them, e.g. xii is freely 
variable. 

132 which is to see by simply inserting the above ansatz: jn = 1, Ui = \* and IA.52I . 

133 The resulting matrix is then A k = Q7 l c( 1 x ^ 1 ) = , , , c ri (f X *£~\ Xk]2 ) = c( Xk ) 
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This gives a total number of I = 2m — l + 2h — 1 = 2(m + h — 1) edges connected 
to every node, since on each longitude circle j are actually 2h nodes: tke m ™ and 
tke nil ~™~ = — the™ 1 ™ , k = l,...h. This also means, that the opposite point of 
Xk on the latitude circle, which is — Xki counts twice in the sum (on both circles). 



The summation on a latitude circle gives: 134 

1 2m- 



E 



Xk - Xi 



2xk 



(A.55) 



The summation on the longitude circle gives (the opposite point of Xk on this 
circle, gives the summand l/2%fc): 



E 

ieA(fe) 



1 



Xk - Xi 



1 

2Xk 



tk_ 

Xk 



E( 

L i^k 



1 



1 



tk t% H~ 



Consequently the complete sum over all edges becomes: 



E 



i 



Xk - Xi 



Xk 

t 2 



E +2 

i^k k 



2*2 



(A.56) 



(A.57) 



To prove corespondence to zeros of spherical harmonics, one subtitutes back to 
the cartesian coordinates, which is given by t 2 = \x\ 2 = fef ■ Then (|A.57|) can 
be expressed as 



E 



Xk ~ Xi 



= Xk( 1 + Z k) 



1 - z k 



E 



Zk 



The r.h.s. of this equation is then equal to the required expression 



l- 



Xk 



l xt 



1 



Zk 



(m + h - l)x* k (l + z k ) 



l + \Xk\ 2 2 
to obey the stationarity conditions l|A.52(l . if for all Zk holds: 



E 



i 



z k 



Zk 



(A.58) 



It is now easy to check, that (|A.58|) is fulfilled for the 
zeros of legendre functions of order m + h, namely P™^(z) (see e.g. 
0, pp. 282), with the methods of orthogonal polynoms presented here (see 
section El eq. I|A.31|I ). 

As a summary of this chapter, I want to state, that all classical stationary 
states of the DIRAC equation are reproduced exactly above. However, it is not 
yet clear, by which principle some grids (e.g. with an odd number of nodes on 
a latitude circle) are suppressed. 

134 This formula can be proved using a general relation for the complex roots of z" = 1, 
which are zj. = e 27rl '=/ n j k = 1, . . . ,n, namely: ^22—1 1— = "2 1 » which again follows from 
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